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5Pi Abstract 
P . 

' We have investigated the properties of in nuclei within the framework of relativistic 

Q>^ ■ mean field. The coupling constants are educed with quark meson coupling model. There 

^—^ , is strong attractive interaction for 0^-nucleus and G"*" can be bind in nuclei. The depth 

^ I of optical potential for G"*" in nuclear matter is estimated. 
> ■ PACS: 21.10.Dr, 26.90.Xh,13.75.Ev. 
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: 1 Introduction 
o ' 

■ Since the extremely narrow positive strangeness pentaquark state Q^{uudds) of mass 1.54 
> . GeV and = 1/2"^ is predicted[l] and with several experimental searches being successfully 

! undertaken[2, 3, 4, 5, 6, 7, 8], pentaquark becomes the hot topic of investigation. The discovery 
^ I of the exotic baryon 0"*" with positive strangeness opens new possibilities of forming exotic 
C ■ 0+ hypernuclei which, like in the case of negative strangeness A, E, S hypernuclei which have 

■ obtained steady progress at the experimental and theoretical levels[10, 11, 12, 13, 14], can 
^ ■ provide information unreachable or complementary to that obtained in elementary reactions. 

Suggestions that 0+ could be bound in nuclei have already been made [15, 16]and the 
selfenergy of the 0^ pentaquark in nuclei is calculated by D. Cabrera eta/. [16], they show 
that the in-medium renormalization of the pion in the two meson cloud of the 0"'" leads to a 
sizable attraction, enough to produce a large number of bound and narrow 0^ states in nuclei. 
However, the spin, isospin and parity of 0^ are still not very determinate, there are many 
discussions about these aspects[l, 15, 17, 18, 19, 20]. In this work, we assume that = 1/2+ 
and / = for 0+ as some suggestions in Refs.[l, 15], thus we can study 0"*" hypernuclei in 
the framework of relativistic mean field (RMF), which has been used to investigate A, S, S 
hypernuclei[21, 22, 23, 24]. 

In the present work, our main aim is to determine the coupling constants a — 0+, uj — 0+ 
for the RMF calculation. Then we attempt to investigate the 0"*" hypernuclei in the framework 
of relativistic mean field, from which we want to know whether the 0"'" baryon could be bound 
in nuclei and the how much the depth of the 0"'"-nucleus potential. 

This paper is organized as follows. The framework of hypernuclei in RMF is introduced 
in Sec. 2. The coupling constants are educed and determined in Sec. 3. The calculations and 
analysis is given in Sec. 4. The last Sec. is the summary. 
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2 Hypernuclei in RMF model 

In the RMF theory, the effective Lagrangian density which includes the nucleon and ©"'"-hyperon 
can be written as 

C — C-^i -\- Ce+ , (1) 



-g'^i^Nl'^plT^'^N + \d^<jd^a - \mla' - \gia' 
1 o . 1 o „ 1 



-gla' - -Q'^'^Q,^ + -m>^a;, - -R^^'^R;^ 



-^®*G+7V^"*G+ + e^e+T^^l^^ll - T')*e+ (3) 

with 

pnu ^ d^'A'-d''A^. (4) 

where the meson fields are denoted by cr, u)^ p^, and their masses by mo-, m^, nip, re- 
spectively; \E'Ar and \E'0+ are the nucleon and G"*" baryon fields with corresponding masses Mn 
and Mq+, respectively; is the electromagnetic field. g^{g^), g^{g^) and gp (gf) are the 
cr — N{a — Q), LU — N{u! — 6) and p — N{p — O) coupling constants , respectively. The isospin 
Pauli matrices are written as r", being the third component of r". For baryon, we 
assume isospin / = 0, r" = 0, there is no coupling to p meson. 



3 How to determine the coupling constants 

According the quark meson coupling model developed in Ref. [25] , the baryon coupling constants 
can be determined. The equation of motion for meson field operators are 

d^d^'a + mla = glqq, (5) 
d.d^cu'' + mluj^ = glq^q, (6) 

d^d'^p^''^ + mlp^''^ = gW^-^q, (7) 

where g^, g^ and g^^ are the quark-meson coupling constants for a, uo and p, respectively. 

The mean fields are defined as the expectation values in the ground state of the nucleus, 
\A): 

{A\a{t,r)\A) = a{r), (8) 
(A|d;^(t,r)|A) = 5(z/,0V(r), (9) 
(^|p'^'"(t,r)|yl) = 5(i/,0)5(a,3)p(r), (10) 



2 



where (t,r) are the coordinates in the rest frame of the nucleus. In the mean field approximation 
the sources are the sum of the sources created by each nucleon - the latter interacting with the 
meson fields. Thus we write 

QQ{t,r) = J2 {Q<li't^r))i, (11) 

i=l,A 

qYq{t,r)^ {qY q{t,r))i, (12) 

i=l,A 

where denotes the matrix element in the nucleon i located at at time t. According to 
the Born-Oppenheimer approximation, the baryon structure is described by n quarks in the 
lowest mode, therefore, in the instantaneous rest frame (IRF) of the nucleon we have 



{q'q'{t',r% = = ns,{u'), (13) 

m 

{qYq'if, r')). = nS{u, 0) ^ 0f '"^(«')0°'"^(«') = riS{u, 0)w,{u'), (14) 



where {t', r') are the the coordinates and q' is the quark field in the IRF. 0°'™ is the complete 
and orthogonal set of eigenfunctions which are defined by Eq.(23) in Ref.[25]. 
Prom the Lorents transformation properties of the fields we get 

{qq{t.r))i = ^((coshC,))"^ / dke^^<^-''^~^^ S{KR,), (15) 

{qi\{t,r)), = ^/ dke^^<^-''^^^W{k,R,), (16) 

{qiq{t. r))i = j^^v, J dke''<^-''^^^W{k, R,), (17) 

with 

S{k,Ri)= j due-'^^^-''^+^^''^'^°^^^'hi{u), (18) 

W{k,Ri) = j diie-^('^^-"^+^^"^/"°^'^^'Vi(ii). (19) 

The corresponding physical quantities such as ^j, ul etc. are defined in Ref.[25]. Finally, 
the mean field expressions for the meson sources take the form 

{A I qq{t,r) \ A) = dke'^^A \ Y,{{cosH,))-'e-'^-''^S{KR,) \ A), (20) 

{A I qYq{t,T) \A)^-^J dke^^^A \ Y^e-'^'^'WiKR,) \ A), (21) 

{A\q^q{t,r)\A)^Q. (22) 
To simplify further, we remark that a matrix element of the form 

{A\Y,(^-'^''^...\A) (23) 

i 

is neghgible unless k is less than, or of the order of, the reciprocal of the nuclear radius. 
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We now define the scalar, baryon and isospin densities of the nucleons in the nucleus by 

,. = {A\^^^iir-B.)\A), (24) 

PB = {A\Y.6(r-Ri)\A), (25) 

i 

P^ = {A\Y.%5{t-R,)\A), (26) 

where we have used Ei = M^cosh^j + V{Ri) (Eq.(58) of Ref.[25]) to eliminate the factor 
(cosh^j)~^ The meson sources take the simple form 

{A\qq{t,r)\A)=nS{T)ps, (27) 

{A\q^%{t,r)\A)^npB, (28) 

{A I Tr^-^q{t,r) I A) = 5(i/,0)5(a,3)p(3). (29) 
We have used the notation 

l]o/2 + m*i?B (1^0-1) 



5(r) = 5(0,r) = j dusr{u) 



no{^o-l)+mlRB/2' 



^0 = + (i?sm*)2, 

m* = mq-gla{r) (30) 

Here, x ^'^^ ^'^^ the parameters for lowest eigenvalues for the quarks and the corresponding 
current quark masses. Rb is the radius of the bag. 

Since their sources are time independent and since they do not propagate, the mean fields 
are also time independent. So, we get the desired equations for (7(r), uj{r) and p(r): 

{-A + ml]a{r) = g,C{r)ps{r), (31) 
{-l^ + ml}uj^{r)=g^PB{r), (32) 
{-A + ml}po^{r)=gpP^{r). (33) 

where the nucleon coupling constants and C(r) are defined by 

g, = nglS{a = 0),g^ = ng^g, = g%C{T) = S{r)/S{a = 0). (34) 

In an approximation where the (T{r) uj{r) and p(r) fields couple only to the u and d quarks[26, 
27], the coupling constants in the strange baryons are obtained as g^ = {nq/n)gi^ and g^ = 
9p = 9% with Uq being the total number of valence u and d quarks in the baryon S. For the 
strange baryon S, the coupling constant a — S can be written as 

g'^=nqglSs{a = 0), (35) 

where 

^^f--' = a„(fi„-l) + m.fl./2- 
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where Rs is the bag radius for the strange baryon. By combining the Eq.(24) and Eq.(25) we 
get 



= %.Ss{a = 0)/S{a = 0) = ^.Ts/b- (37) 



According Eq.(24) 



9a = -irQa : 9u, = ^^o; : 9p ^ 9p : (38) 



where , g^ , g^ are the couphng constants for nucleon. So, we can easily get the relations 

„5_"'9_Arp „S_''^q^N S_ N foQ\ 

9a - yfl'a i s/B, 9u> - -^9uj ^9p - 9p ■ (3yj 

4 Calculations and Analysis 

In this section we will show the our numerical results using relativistic mean field theory model. 
In the mean-field and the no-sea approximations [28] the contributions of anti-(quasi) particles 
and quantum fluctuations of mesons flelds can be thus neglected. The contribution of the tensor 
coupling is neglected for simple, although it is considered in some Refs.[22, 23]. The coupling 
constants of nucleon to mesons employed are NL-SH[30], which are list in table 1. The couphng 
constants of ©"•" to mesons are obtained from the relations Eq.(39). For ©+, ng=4, isospin 7 = 
means no coupling with p meson. Thus, we have the relations 

9t-\9^^e/B,g^^=\g!:. (40) 

In practice the value for Ts/s ~ 1 for all hyperons even though Rb 7^ -R5[29], thus in our 
calculation, we use Tqjb — 1- 



4.1 The results for 0+ hypernuclei in RMF 

In flgure 1, wc show the single-particle energy spectrum such as lsl/2, lp3/2, lpl/2, ld5/2... 
of O"^ in ^Li,^^C, ^^O, "^^Ca and ^°^Pb, which is denoted by 1, 2a, 3, 4a, 5 respectively. 2b is 
the energy spectrum for @ C when the coupling constant change a little to g^ = ^g^ x 0.95. 
2c, 2d and 4b, 4c are the results for q C and Ca at the depth of potential 60 MeV and 120 
MeV respectively of rcf[16]. 

In light nuclei like ^Li (1), ^^C(2a) and ^^0(3) we flnd several bound states. The separation 
of the two deepest bound states (lsl/2 and lp3/2) is about 24 MeV for ^Li, about 60 MeV 
for ^^C and around 30 MeV for ^^O. For medium and heavy nuclei, as in '^'^Ca,^°^Pb shown 
in the figure, we find more bound states and the energy separation between lsl/2 and lp3/2 
is about 18 MeV for ^°Ca and about 4 MeV for ^°^Pb. There is a trend that the separation 
between lsl/2 and lp3/2 becomes smaller and smaller and the bound state becomes more and 
more with the increasing of the nucleon number. It is interesting that the width of G"*" is very 
narrow with upper limit as small as 15 MeV predicted by the chiral soliton model in Ref. [1], 
and 9 MeV [3 7]. Studies based on K'^N scattering suggest that the width should be smaller than 
5 MeV[34, 35, 36]. Cahn and triUing have extracted r(e) = 0.9 ± 0.3 MeV from an analysis 
of Xenon bubble chamber [38]. The narrow upper limit width of O"*" is often smaller than the 
separation of the two deepest bound states (lsl/2 and lp3/2) for light and medium nuclei, 
which would make a clear case for the experimental observation of these states. 
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The splitting of spin-orbit for ©+ hypernuclei can be get from our calculation with RMF. 
Such as, from figure 1, we can find that the splitting between lp3/2 and lpl/2 is about 5 MeV 
for ^Li (1) and ^^O , 10 MeV for ^^C, The splitting is smaller for medium nucleus '^''Ca (about 
2 MeV) and there is no obvious splitting for heavy nucleus ^°^Pb. If the prediction based on 
K^N scattering suggest that the width should be smaller than 5 MeV[34, 35, 36] is true, the 
the splitting between lp3/2 and lpl/2 is also bigger than the narrow upper limit width of G"*", 
which means that we may observe splitting of spin-orbit for light hypernuclei. There is a 
very interesting phenomenon that both separation between lsl/2 and lp3/2 and splitting of 
spin-orbit between lp3/2 and lpl/2 are much larger than the others. 

For Comparison, we also yield the the single-particle energy spectrum of 0+ from Ref[16] 
by solving Schrodinger equation with two potentials. It is a very simple model and can not 
give the splitting of spin-orbit. It estimates that the binding energy(ls) of ^^C is in the region 
34.0 ~ 87.3 MeV and 42.6 ~ 98.2 MeV for ^"^Ca. In our calculation, the binding energy(ls) 
of ^^C is about 116 MeV and 81 MeV for ^°Ca, and the energy separation between two bound 
states is about 18 MeV for '^°Ca and 60 MeV for ^^C. Our results are very close to the the 
results estimated by solving Schrodinger equation in ref[16]. 

In table 1 we list the calculated binding energy per baryon, —E/A, r.m.s. charge radius, rch, 
and r.m.s. radii of the O"*" baryon and the neutron and proton distributions (re , r„ and r^). 
For comparision, we also give these quantities without 0+. The differences in vulues for finite 
nuclei and hypernuclei listed in table 1 reflects the effects of the 0+. The appcarcncc of 0"*" 
in nuclei leads to larger binding energy per baryon, —E/A. From table 1, it is found that the 
values of r^, Vp and rch for eBe, @ C are obvious less than those for the corresponding normal 
nuclei. This is the shrinkage effect which has been found in lighter A hypernuclei [24, 31]. In 
our calculation for ©Be, g C, we also flnd that there is shrinkage effect and gluelike role for 0"*". 
But this phenomenon is not found for @Li and it is not obvious for qO, q Ca and Q^^Pb. 

The key point for our calculation is how to determine the coupling constants of 0^ to 
mesons as accurate as possible. In the calculation, we find that the binding energy is very 
sensitive to the coupling constants. In Refs[21], they determine the the coupling constants of 
hyperon to mesons by adjusting / around 2/3 (according Eq.(39),5'J/(?^ ~ 2/3 ) to flt the 
experimental data because of the effect of medium in nuclei. Where Y is stands the hyperon. To 
consider the effect of medium to the coupling constant 5^®, we list the result of ^^C in figure 1, 
denoted with 2b, when the couphng constant changes to g'^xO.QS. From the figure, we see 
the binding energy of lsl/2 have a large change from 116 to 60 MeV, but there is still several 
bound states and the energy separation energy between lsl/2 and lp3/2 is around 32 MeV. 

We also the list results of the ground state for ^^C when the coupling constant (7® changes 
from 13.925 to 11.725 (c/®x0.84) in table 2. From it, we can see that the binding energy of 
lsl/2 is sensitive to the coupling constant, if g^ < 11.725(13.925 x 0.84), there are no bound 
states for 0+ in ^^C. In Ref.[21], g^/g^ is 0.621 (0.667 x 0.93) and 0.619 (0.667 x 0.928) for 
A — o" and S — o", respectively. Thus, from the calculation of before[21],we can conclude that the 
effect of medium is not very large, the considered medium effective coupling constant is about 
0.93 x g^ for hyperon. If we let gf® = 0.93 x 13.925=12.95, from table 1, we can see that the 
binding energy of lsl/2 for ^^C is between 52 and 40 MeV, which is agree with the estimation 
34 ~ 87 MeV in the Ref[16]. 

In table 2, the binding energy per baryon, -E/A (in MeV), r.m.s. charge radius Vch (in fm), 
r.m.s. radii of O"*", neutron and proton, tq , r„ and Vp (in fm) are also shown respectively. It 
is seen in table 2 that the radii is increasing and the binding energy per baryon is decreasing 
with the decreasing of the coupling constant. When gf < 12.325, re < rch, which means that 
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the ©+ is out of the nucleus ^^C. It can only form ©+ atom. 



Table 1: The parametrization of the nucleonic sector (adopted from Ref.[30]). The masses are given in MeV 
and the couphng §2 in fm~^. 





939.0 


9^ 


10.444 


93 


-15.8337 


rria 


526.059 


9!i 


12.945 








783.0 




4.383 






rrip 


763.0 


92 


-6.9099 







Table 2: Binding energy per baryon, —E/A (in MeV), binding energy of lsl/2 for in nuclei, 
—B (in MeV), r.m.s. charge radius Tch (in fm), r.m.s. radii of hyperon, neutron and proton,re 
, r„ and (in fm), respectively. The configuration of 0+ is lsl/2 for all hypernclei 





rp(fm) 


r„(fm) 


r0(fm) 


rc/i(fm) 


-E/A{McN) 


-S(MeV) 


«Li 


2.37 


2.32 




2.51 


5.67 




^Li 


2.36 


2.36 


1.38 


2.50 


6.21 


52.13 


8Be 


2.34 


2.30 




2.48 


5.42 




|Be 


2.21 


2.19 


1.35 


2.36 


11.43 


67.72 


12(3 


2.32 


2.30 




2.46 


7.47 




13p 


2.17 


2.15 


1.19 


2.32 


14.41 


116.24 


160 


2.58 


2.55 




2.70 


8.04 




e ^ 


2.58 


2.54 


1.46 


2.70 


8.60 


83.66 


^"Ca 


3.36 


3.31 




3.46 


8.52 




11 r-., 

(-) 


3.35 


3.30 


1.86 


3.15 


10.18 


81.35 


208pb 


5.45 


5.71 




5.51 


7.90 




|°9Pb 


5.44 


7.70 


3.87 


2.49 


8.19 


69.54 



4.2 potential depth in nuclear matter 

In nuclear matter the potential depth is written as [24] 



(41) 



where o"^* and uj^ are the values of a and uj fields at the saturation nuclear density, respectively. 
Use the relations cr^* = (Mj^ — MN)/g^ and loq'^ — Qo/mlj at the saturation nuclear density, 
we can get [32] 



.Ml 



.9^ 



-.N 



1© 



AT 1 \ ^0- , / 9uj \2 9uj ^ 



(42) 



'Mm ' gi; m^' gi 

where M^ is the effective mass of nucleon and Qq is the saturation nuclear density. From Eq.(42), 
we can see that the depth of potential U is determined by ^p^, % and if we let % = 4/3 



Table 3: Binding energy per baryon, —E/A (in MeV), binding energy of lsl/2 for 0"*" in nuclei, 
—B (in MeV), r.m.s. charge radius Vch (in fm), r.m.s. radii of hyperon, neutron and proton,re 
, r„ and Vp (in fm), respectively. The configuration of is lsl/2 for @C in different g^. 



9^. 


rp(fm) 


r„(fm) 


re(fm) 




-E/A{MeV) 


-S(MeV) 


13.925 


2.17 


2.15 


1.19 


2.32 


14.24 


116.24 


13.725 


2.18 


2.16 


1.22 


2.33 


13.27 


100.00 


13.525 


2.19 


2.17 


1.28 


2.34 


12.33 


83.70 


13.325 


2.21 


2.19 


1.35 


2.36 


11.43 


67.72 


13.125 


2.24 


2.21 


1.46 


2.38 


10.60 


52.92 


12.925 


2.26 


2.24 


1.59 


2.41 


9.85 


40.17 


12.725 


2.29 


2.27 


1.76 


2.44 


9.19 


29.67 


12.525 


2.32 


2.29 


1.95 


2.46 


8.60 


21.20 


12.325 


2.34 


2.31 


2.20 


2.48 


8.10 


14.43 


12.125 


2.35 


2.32 


2.50 


2.49 


7.68 


9.08 


11.925 


2.35 


2.32 


2.90 


2.49 


7.34 


4.98 


11.725 


2.35 


2.32 


3.43 


2.49 


7.10 


2.01 



as in Eq.(40). Usually, the effective mass -g^ and saturation nuclear density are considered 

as "experimental input". In the RMF calculation, ~ 0.60, Qq ~ 0.15fm^^[33]. Here, we 
let them have some uncertain and can change in a little region. According the calculation and 
analysis of section 3, when < 11.725(13.925 x 0.84), there are no bound states for ©+ in 
^^C. thus we let it as a approximation that there is no bound states for 0+ in nuclear matter 
when 5® < 11.725(13.925 x 0.84), namely the depth of potential C/ > . 

M* 

In figure 2, we show the relation of the potential depth changing with the effective mass 
for qq = 0A5,gQ = 0.16 and qq = 0.17fm~^ in a,b and c respectively. From figure 2(a), it is easily 
seen that U=0 corresponds to effective mass = 0.60 when qq = 0.15fm~^. With the same 
method we can determine = 0.575 for qq = 0.16fm~^ and = 0.546 for qq = 0.17fm~^ 

from b and c. The results of effective mass is reasonable, because they are around 0.6, 
which agrees with the fit data in RMF [33]. 

M* 

Since the effective mass and saturation nuclear density ^-re determined, we can get the 
potential depth at saturation nuclear density easily. From figure 2 (a), we can see U—-S1 MeV at 

= 13.925 and f/=-46 MeV at gf = 13.925x0.93 when ^ = 0.60,^o = 0.15fm-=^. Considered 
the effect of medium, if the (7—0+ coupling constant is in the region 13.925 < gf < 13.925x0.93, 
the potential depth is in the region —81 < U < —46 MeV. From figure 2(b) and (c), we also 
get -84: < U < -47 MeV for = 0.16fm-3 and -92 < U < -52 MeV for = 0.17fm-3 
when 13.925 < gf < 13.925 x 0.93 in the same way. From the calculation, it is found that 
the uncertain of the saturation nuclear density can bring some uncertain for the potential 
depth in several MeV. Because there is still no experimental data for 0"*" in nuclear matter, 
we can only estimate the potential depth in theory, the uncertain for the potential depth in 
several MeV is allowed. In fact, the potential depth estimated by us agrees well with that of 
in Ref.[16]-120 < U < -60 MeV. 
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5 Summary 



We have educed the coupUng constants for strange many quarks system with quark meson 
couphng model. We get the relation gf = ^Q^Tq/b , = ^g^ for 0+ coupling with a and 
LU meson. With the coupling constants for a — and uj — educed by us, the calculations 
for 0+ in nuclei are carried out in the framework of RMF. Prom the single-particle levels of 
0"^ in nuclei, it is found that there are several bound states for light nuclei, such as @Li, and 
more bound states for heavier nuclei. The binding energies of lsl/2 are usually as high as 50 
MeV to 120 MeV, which means that there are strong attractive interaction for ©^-nucleus. This 
prediction of RMF is agree with the prediction of Refs.[15, 16]. The strong attractive interaction 
for ©^-nucleus induces stronger per nucleon binding energy than normal nuclei, which can be 
seen obviously in table 2. Because of the strong attractive interaction for ©^-nucleus, the 
shrinking effect is found for light nuclei qBc and @ C, in medium or heavy nuclei there is no 
obvious shrinking effect. The shrinking effect is first found in the light A-hypernucleus for 
lithium hypernuclei[30], but there is no shrinking effect for qLI in RMF prediction. 

The depth of optical potential for ©+ in nuclear matter is investigated, the region is given. 
If we assume the saturation nuclear density = 0.15 0.16 and 0.17 fm~^ with the couphng 
constant 13.925 < gf < 13.925 x 0.93, the depth of potential -81 < U < -46 MeV, -84 < 
U < -47 MeV and -92 <U < -52 MeV, respectively. 

In this work we assume that = 1/2"*" and 1 = for ©'^, and when we induce the coupling 
constants in section 3 with QMC model, we do an approximation where the a{r) iL!{r) and p(r) 
fields couple only to the u and d quarks. These assume and approximation may be affect our 
results. 
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Figure 1: Single-particle levels for in ^Li, ^^C, ^^O, "^^Ca and ^osp]^ denoted with 
l,2a,2b,2c,2d,3,4a,4b,4c and 5 respectively. 1, 2a, 3, 4a, 5 are our results with RMF at 
gf = 13.925 and 2b is our result at gf = 13.925 x 0.95. 2c,2d and 4b,4c are the results 
from Ref. [16] in ^^C and ^°Ca at the depth of potential 60 MeV and 120 MeV respectively. 
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Figure 2: The depth of potential for 0+ in nuclear matter versus effective mass M* /M for 
fixed values of U at the saturation nuclear density. a,b and c is for = 0.15, 0.16, 0.17fm~^ 
respectively. The dash dot line is for gf = 13.925 x 0.84, dash fine is for gf = 13.925 x 0.93, 
and solid fine is for gf — 13.925. 
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